Exact Relativistic Ideal Hydro dynamical Solutions 
in (1+3)D with Longitudinal and Transverse Flows 
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A new method for solving relativistic ideal hydrodynamics in (1+3)D is developed. Longitudinal 
and transverse radial flows are explicitly embedded and the hydrodynamic equations are reduced 
to a single equation for the transverse velocity field only, which is much more tractable. As an 
application we use the method to find analytically all possible solutions with power dependence on 
proper time and transverse radius. Possible application to the Relativistic Heavy Ion Collisions and 
possible generalizations of the method are discussed. 
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I. INTRODUCTION 



Relativistic hydrodynamics has wide applications in 
a variety of physical phenomena, ranging from the 
largest scales such as in cosmology and astrophysics [1] 
to the smallest scales such as in the relativistic nuclear 
collisions 0, For an introduction on the general for- 
malism, see e.g. d, m, 

Recently there has been a remarkably successful ap- 
plication of Relativistic Ideal Hydrodynamics (RIHD) to 
the description of the space-time evolution of the hot 
dense QCD matter created in the Relativistic Heavy Ion 
Collider (RHIC) experiments. In the collisions of two 
relativistically moving heavy nuclei, a lot of energy is 
deposited in a small volume which soon creates an equi- 
librated system of high energy density with special ini- 
tial geometry: extremely thin in the beam direction z 
while in the transverse plane x — y it is of the size of 
the nuclei. The space-time evolution at RHIC is charac- 
terized by fast longitudinal expansion (longitudinal flow) 
and relatively slowly developing transverse expansion (ra- 
dial and elliptic flow). In non-central collisions the cre- 
ated matter on the transverse plane x — y is initially 
anisotropic: such initial spatial anisotropy leads to differ- 
ent pressure gradients and thus different accelerations of 
the flow along different azimuthal directions. The result- 
ing anisotropic transverse flow velocity eventually trans- 
lates into the anisotropic azimuthal distribution of the 
final particle yield which is represented by the experi- 
mental observable called elliptic flow V2 — one of the 
milestone measurements at RHIC [8]. The RIHD model 
calculations [1, [1(1 HH (and more recently its extension 
to include viscous corrections fl^b performed with re- 
alistic initial conditions and Equation of State (E.o.S) 
for RHIC, are able to reproduce the elliptic flow data 
at low-to-intermediatc transverse momenta for almost all 
particle species and for various centralities, beam ener- 
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gies and colliding nuclei. These achievements of RIHD 
have been the basis for the RHIC discovery that the 
matter b eing created is a strongly-coupled nearly-perfect 

fluidd M HE M 

with extremel y s hort dissipative 
length. It has been suggested [IE [3 III, III [III] that the 
microscopic origin could be due to the strong scattering 
via Lorentz force between the electric and magnetic de- 
grees of freedom coexisting in the created matter, with 
the magnetic ones ultimately connected with the mecha- 
nism of QCD confinement transition. 

The great success of RIHD at RHIC has also inspired 
considerable interest in the formal aspects of relativistic 
hydrodynamics, particularly in analytical solutions of the 
RIHD equations with an emphasis on possible application 
to RHIC, see e.g. [H, H M, M, M, ffli 11 HI . The idea 
to use exact simple RIHD solutions to describe the multi- 
particle production in high energy collisions dated back 
to Landau and Khalatnikov [30| . An important solution 
came from Hwa and Bjorken's works [U [Hj], i-e. the ra- 
pidity boost invariant (1+1)D solution which is widely 
used to describe the longitudinal expansion at RHIC. 
Many of the above mentioned recent works [23ll24l. [25l[26j 
concentrate on finding (1+1)D solutions that give an al- 
ternative description of the longitudinal expansion and 
a more realistic (non-boost-invariant) multiplicity distri- 
bution over rapidity. 

Despite the progress in solving RIHD in (1+1)D, it is 
quite difficult to solve them in higher dimensions. To 
develop methods and find solutions in a realistic (I+3)D 
setting with potential application for RHIC remains an 
attracting but demanding task. In this work, we will 
develop a new method to find solutions in (I+3)D with 
both longitudinal and transverse flows. In Section III, 
we will show how the method can reduce the hydrody- 
namics equations to a single constraint equation for the 
transverse velocity field only. Using the derived equation, 
we will find all solutions with power-law dependence on 
proper time and transverse radius in Section IV. The 
physical relevance of our results to RHIC and possible 
generalizations will be discussed in Section V. We will 
also include a brief introduction of RIHD in Section II 
and an illustration of the method in (f+l)D in the Ap- 
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pendices B and C. 



II. REVIEW OF RELATIVISTIC IDEAL 
HYDRODYNAMICS IN (1+3)D 

The hydrodynamics equations in general are simply 
the conservation laws for energy and momenta, i.e. 



T mn , „ = 



(1) 



with to, n running over 3+1 space-time indices. Fol- 
lowing usual convention (in e.g. 0, [1]), the subscript 
"; n" denotes the covariant derivative D n while a sub- 
script ", n" is for ordinary derivative d n . Below we will 
introduce curved coordinates in order to simplify the hy- 
drodynamics equations. Therefore, we use the general 
form for the hydrodynamics equations which involves co- 
variant derivatives, see e.g. [HIH ■ Throughout this paper 
we discuss only hydrodynamics without any conserved 
charge, leaving the situation with conserved currents for 
further investigation. 

For relativistic ideal hydrodynamics, the stress tensor 
is given by 



(e + p) u m u r ' 



pg 



(2) 



with e,p the energy density and pressure defined in the 
flowing matter's local rest frame (L.R.F) which by defi- 
nition are Lorentz scalars. The flow field u m (x) is con- 
strained by u m ■ u m = 1- In the usual {t, x) coordinates 
one can express u m (x) as 7(1, v) with 7 = — W and 

v = dx/dt. 

We further need to specify an equation of state (E.o.S) 
relating the energy density e and the pressure p of the 
underlying fluid. Here we employ a simple, linear E.o.S, 
which is typic ally used in analytic studies of RIHD (23l . 

El, HE M M Hi: 



p = v(e+p) 



(3) 



The above means e 



= ^—^-p implying a speed of sound 
and, in order to assure c s < 1, we 



de — Y 1- 

require < v < 1/2. We note, that by adding a constant 
to the above E.o.S, the resulting velocity field remains 
unchanged. 

The hydrodynamics equations together with the E.o.S 
thus form a complete set of 5 equations for the 5 field 
variables: e(x),p(x) and the three independent compo- 
nents of u m (x) . 



A. Hydro Equations in Curved Coordinates 

When formulating hydrodynamics for application to 
e.g. the relativistic heavy ion collisions, it is often useful 
to use alternative coordinates systems which are curved. 
For our purpose of studying the (1+3)D solutions with 



longitudinal and transverse flow, we will use a coordinate 
system of (r, rj, p, (f>): i.e. the proper time, the spatial 
(longitudinal) rapidity, the transverse radius and the az- 
imuthal angle. They are related to the usual (i, x, y, z) 
in the following way: 



•\A 2 — z 2 , 77 = -In 



1. t + z 



2 t-z 1 

o / 1 1 x + y ■ i 

p = \Jx l + y z , (/) = —In : 

li x — y ■ 1 



(4) 



and inversely 



t = t cosh 77 , z = t sinh 77 , 
x = p cos 4> , y = p sin <f> 



(5) 



The velocity field u m in these coordinates is related to 
= 7(1, v) in flat coordinates (f, x) via 

7 

u T = 7 (cosh rj — v z sinh 77) , u v — —(v z cosh 77 — sinh 7/) , 

T 

. 7 

u p = j{v x cos 4> + v y sin cp) , v? = — (v y cos 4> — v x sin <f>). 



(G) 



The metric tensor associated with the (t, 77, p, (j)) coor- 
dinates is 



9mn = Diag(l, -t 2 , -1, -p 2 ) 
g mn = 13^(1,-1,-1,-1) 



(7) 



For the covariant derivatives we will need the Affine 



m , n ~r Qkn , m 



connections T 3 mn = g jk T kmn = g jk \{gh 
gmn . k) ■ In our case, the non- vanishing connections are: 

1 



F T = r T v = r v 

VV ' nr T7J 

pp _ _ r _ r 4> _ 1 



(8) 



We also give the explicit form of covariant derivatives in 
the present coordinates for an arbitrary contra-variant- 
vector A k (i.e. with upper indices k — r,rj, p, <j>) : 



A k , T 


= A k , T - 


(- T k i A i 




= 


i-r*^ 


Ak 


= A k . p - 


VT k pi A* 




= A k ^- 


r-T^A 1 



t 5 k A' 1 + -S k A T 

T 1 

-si A* 



pS k A+ + -6i AP 



(9) 



Inserting the above Eqs. (|2l6|7|) into the general hydrody- 
namics equations ([TJ) and making use of Eq. (|8l9p . one 
obtains the hydrodynamics equations explicitly in the 
curved coordinates (r, 77, p, 4>). 
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B. Some Known Simple Exact Solutions 

We now recall some known simple exact solutions that 
are pertinent for our approach. 

One of the most famous examples is the so-called Hwa- 
Bjorken solution [3l], H2] which is essentially the Hubble 
expansion in (1+1)D. The pressure and velocity fields of 
this solution are given by 



constant 

PBj. = t1/{1 _ v) 

u Bj . = (1,0,0,0) 



(10) 



It is more transparent to look at the components of v in 
flat coordinates, which are simply v z = tanh?7 = j , v x = 



Vy = 0. 



A generalization of the Hwa-Bjorken solution to radial 
Hubble flow in (1+3)D (with further generalization to 
(l+d)D ) is straightforward. The pressure and velocity 
fields are 



constant 



PHu. 



(t 2 



n 2 ) 3(1-1/) 



cosh 77 



uhu. = ^( CO g h '°' ~'°)< 7= — j — 

^ T \Jl — (p/t) 2 cosh 2 n 

(11) 

In flat coordinates the velocity fields are given in simple 
form: v = x/t = (x/t,y/t, z/t). 



III. THE NEW REDUCTION METHOD 

In this section, we use a new reduction method to 
find solutions for (1+3)D RIHD equations. The general 
idea is to first embed known solutions in lower dimen- 
sions which automatically solve 2 out of the total of 4- 
component hydrodynamics equations, and then reduce 
the remaining 2 equations into a single equation for the 
velocity field only. As usual, one starts with a certain 
ansatz for the flow velocity field: in our case we will 
use an ansatz with built-in longitudinal and transverse 
radial flow, aiming at possible application for RHIC. It 
would be even more interesting to include transverse el- 
liptic flow which requires a suitable curved coordinates 
(like certain hyperbolic coordinates) other than the one 
used here. However generally in those cases, more Affinc 
connections are non- vanishing, which makes the reduc- 
tion method discussed below much more involved: we 
will leave this for future investigation. 



A. Including Longitudinal and Transverse Flow 

We first embed the boost-invariant longitudinal flow as 
many numerical hydrodynamics calculations do, which is 
a suitable approach for RHIC related phenomenology. To 
do that, we simply set v z = z/t = tanhr/, i.e. u n = 0. 



Next we include the transverse radial flow which is 
isotropic in the transverse plane. Radial flow is substan- 
tial and important at RHIC. To do so, we introduce the 
radial flow field v p and set the flat-coordinate transverse 
flow fields to be v x — v p cos (ft and v y — v p sin (ft, which 
implies for the curved coordinates u p = jv p and = 0. 
We note that this ansatz goes beyond a simple change 
to cylindrical coordinates, since we require that = 
which considerably simplifies the hydrodynamics equa- 
tions. 

To summarize, in order to describe a situation with 
both longitudinal flow and transverse radial flow we have 
made the following ansatz for the flow fields u m in the 
coordinates (t, n, p, (ft): 

u m = j(l, 0, v p , 0) (12) 
v p = v p cosh n , 7 = 1/ yjl-v 2 

Note that we need to require v p < 1. 



B. The Equation for Transverse Velocity 

With the flow fields given in (|12p. we can now explicitly 
express the stress tensor components. The non-vanishing 
ones are given below: 

T TT = f{e+p)-p=C-~l)p (13) 
v 

T pp = fv 2 p {e+ V )+p=C-^ L + l)p (14) 



T T <> = fv p {e + p) = ^p 



rp4><j> 



V 
P_ 



(15) 
(16) 



For the second equalities in each of the first three lines 
we have used the E.o.S ([3j> to substitute e+phy p/v. 

With the above expressions and using (jHJ ([9]), the hy- 
drodynamics equations (fT]) then become 



T TT T + + " + T t p p + = 0(17) 

T T p 



(18) 



T pX . A = T pp , p + + T Tp . T + = 0(19) 



rprp 



P 



(20) 



The two equations involving derivatives over n and (ft are 
trivially solved by setting p(x) — p(r, p) (and the same 
for energy density e(r,p) due to the E.o.S) and accord- 
ingly v p (x) — v p (r,p). We note, that the simple form 
of Eqs. (|18l20[) . are a direct consequence of the vanishing 
components u v — = in the flow field ansatz, Eq.(fT2]). 
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Finally we introduce a combined field variable K, de- 
fined as 

T TT +p 7 2 p vpr 
K, = — — — = > p = -^-K (21) 



We then substitute the pressure p in the equations (|17|) 
and (Til?|) and obtain two equations for the fields K. and 
v p , which can be expressed as 

D a -fC, T +D b -fC, p = Dx-K (22) 
D b -K. lT + D c -IC l p=D 2 -K (23) 

The coefficients D a , D b , D c , Di,D 2 are given by: 

D a = (1 - v) + wij 
-Dfa = w p 

D c = V + (1 - 

Di = -IvVpVp , T - v(l - v 2 )/t -v PtP 

D 2 = -2(1 - v)v p Vp iP + 1/(1 - v 2 )/p -v p , T 



(24) 



From [|22"l)(l23l we obtain 



IC,t n r , DcDi - D b D2 _ , 

= ( ln ^),r = -phf; — Fir- = jr M 



/C 



/c 



a; 



^ = (ln/C), p = 



D a D c - D 



2 =^[ T >P] 



(25) 



with the functions given by 

1 



[(l-v)Vp-v]v PjP + 
\{l-2v 2 ) + 2v{v-l)v 2 p ]vpVp, 
-u[l - v%v + (1 - V )v 2 p ]/t - vv p [l - vl]/p\vHV) 

S[vp(r,p)) 



v(l-u)(l-vj)* 



[{vv 2 p + (v-l)}vp >T + 
[(-1 + 4i/ - 2v 2 ) + 2v(y - l)v 2 p ]v p v p , p 



+vv p [l - v 2 ]/r + u[l - v 2 }[{l -v) + vv 2 ]/ p \ (27) 



In (|25|) . the function ln/C depends (via v p ) on two vari- 
ables r and p, and we have two equations for the two 
first order derivatives 91 " ^ and P K . For In /C as a sin- 
gle function of two variables r, p, the two equations can 
be consistent only if the following constraint on second 
order derivatives are satisfied d a ^g p K ' — d g p J c ]J C , i- e - 

d d 

—F[vp( T ,p)]-—g[vp(T,p)} = (28) 



Thus we only need to solve the above single equation 
for the velocity field v p (t, p). Since J- ', Q already involve 
the first derivatives of v p jT and v PtP , the reduced ve- 
locity equation [28] is a second-order partial differential 
equation for the velocity field. As a minor caveat, the 
method applies to the case/region in which ln/C is at 
least second-order differentiable. This reduction method 
can be demonstrated in the more explicit case of (1+1)D 
hydrodynamics, see Appendices B and C. 

Given the above constraints, we can then solve from 
the matter field S directly 



K, = /Co 



_/; o dr'nr',p]+S p P dp'g[T,p']\ 



(29) 



with /Co being the value at arbitrary reference point 
to , Po- 

Finally let us summarize our approach: after including 
into the flow field ansatz the physically desired longitu- 
dinal and transverse flows, we have reduced the hydro- 
dynamic equations into a single equation (|28[) involving 
CWLY~the transverse velocity field v p , and any solution 
to this equation automatically leads to the pressure field 
which together with the velocity field forms a solution to 
the original hydrodynamics equations: 



constant x 



pr_ 

- 9 



x e 



[j T dr'r[r',p] + J" dp'g[T,p'}] 



(30) 



C. Examination of the Method 

We now examine the correctness of the reduced equa- 
tion (|28|) and the solution (|30|) , using the two known sim- 
ple analytic solutions (fT0|) and (fTTjl as both of them are 
certain special cases of our embedding with longitudinal 
and transverse radial flows. 

For the 1-D Bjorken expansion, we have v p Bj. = 
which leads to 



1 



It 



1 

P 



(31) 



One can easily verify that the above Tsj. , Q Bj. satisfy 
the reduced equation (|28[) . Furthermore by inserting 
Tsj. j Qbj. into the solution (|30|) one finds exactly the 
pressure in (TiTj| . 

For the 3-D Hubble expansion, we have v p Hu. = p/ T 
which leads to 



Fhu — — 



3 v — 5/2 2r 



1 — V T 2 — p 2 



1 v-b/2 -2p 

yHu. — 1 : 5 , 

p 1 — v r z — ycr 



(32) 



Again it can easily shown that the above Thu. > Ghu. sat- 
isfy the reduced equation (|28p . Furthermore by inserting 
Fhu.tGhu. into the solution ([30"]) one finds exactly the 
pressure in (|11[) . 
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IV. APPLICATION OF THE METHOD 

As an example for an application of the embedding- 
reduction method in the previous section, we show how 
to find all possible solutions with the following ansatz for 
the radial velocity field: 



A ■ r B ■ P c 



(33) 



with A,B,C arbitrary real numbers. We note that the 
two known exact solutions we mentioned are special cases 
of the above form: the 1-D Bjorken expansion corre- 
sponds to A — while the 3-D Hubble expansion cor- 
responds to A = 1,B = — 1, C = 1. The velocity field 
}, when put into (|2"6"))(|27p . gives the following 



1-25 B + (1 -AB)v + 2Bv 2 



1 



C+(1-C)v 
C(l - 2u) 



1/(1 - v) 
5(1 - 2v) 



p[l - v 2 ' 

v„ 



+ 



r[l 
1 



1/(1-!/) T[l 



1/(1-1/) p[l-«2]2 

-2C -v/(l- v ) 



(34) 



S- 1 

~ 1 - i/ 
B(2i/- 



2CV - j/ - C 



P v(l-v) 
C{2v-l) 1 



-[!-«*] 



1/(1 -f) P[l-^] 2 



i/(l - i/) 



,1212 



(35) 



We have used w P , P = v p ■ C/p and v p tT — v p ■ B/t. As a 
check of the above result, one can verify that by setting 
A = they reduce to ([31]) while by setting A = 1, B = 
— 1, C = 1 they reduce to (|3"2"]) . as they should. 

By inserting (|3"4"|) ([33)1 into equation (|28p , one obtains a 
rather complicated constraint equation for the constants 
A B, C. However after a lengthy calculations, all pos- 
sible combinations of A,B,C solving the equation can 
actually be exhausted. Leaving the detailed (and techni- 
cal) derivations to the Appendix A, we only list the final 
results here: 

• Solution-I: A = with < v < \ (1-D Bjorken) 
- see (El; 



• Solution-II: A = 1,B = — 1,C = 1 with 
< v < \ (3-D Hubble) — see (fTTj): 



• Solution-Ill: A = 1, B = 1, C = -1 with < z/ < 
^ — the solutions are 



x t 2 — z 2 



y t 2 - z 2 



x + 9 i 9 ' V 4- 9 i 9 ' z 

constant 

( T/9 )l/(l-^( p 2 _ r 2)(l-3i/)/(2i/-2^) 



(36) 



• Solution-IV: A = 1, B = 1/3, C = -1/3 with 
z/ = 1/4 — the solutions are 



_ X ,t 2 — Z 2 s2/3 _y ,t 2 — z 2 , 2/3 _ ' 



p = constant x 



( p 2/3 _ r 2/3)2/3 

(H^TS 



(37) 



• Solution-V: A = -1, B = -1, C = 1 with i/ = 1/2 
— the solutions are 



—x —y z 

1 Vy = , V z = — 



p = constant x (r 2 — p 2 ) 



(38) 



It can be verified that these solutions obtained by the 
method introduced here are indeed solutions of the origi- 
nal hydrodynamics equations |T]) . One should notice the 
different applicable kinematic regions in each of the above 
solutions which comes from the constraint that the flow 
velocity shall be less than the speed of light. For our so- 
lutions, the constraint is z < t & x 2 + y 2 + z 2 < t for 
Solution-I,II,V , and for Solution-Ill, IV the constraint is 
z < t & yj x 2 + y 2 + z 2 > t. For a detailed discussion 
about solutions in different regions with respect to kine- 
matic light cone, see e.g. appendices of [23]. 

We notice that all the solutions (except the trivial 
Solution-I with A = 0) satisfy two features (1) B = —C 
and (2) \A\ = 1. The first feature may be due to di- 
mensional reasons. The second feature, \A\ = 1, may be 
heuristically understood in the following way. We first 
consider the case B = —C < 0, i.e. v p = A{p/r) c 
with C > 0: in this case the solution exists in the re- 
gion p < t ■ |A| _1 / C ', and in particular p = for r = 0. 
Thus for any t > 0, the flow front which travels with the 
speed of light, \v p \ = 1, is located at p = r and hence, 
\A\ = 1. Next we consider the case B = —C > 0, i.e. 
Dp = A(r/p) B with B > 0: in this case the solution ex- 
ists in the region p > t-A 1 / 3 with A > 0, separated from 
an empty region by the boundary at p — r • A X I B . At this 
boundary, the flow velocity approaches the speed of light 
v p — ► 1 which enforces the matter density to drop to zero 
in order to avoid an infinite T mn (due to the 7-factor in 
Eq.@). We imagine that at time r = the matter fills 
the whole space and then starts to flow outward, thus the 
boundary also moves outward from the origin with the 
speed of light v p — 1. This again implies the boundary 
should lie at p = t requiring A = 1. 

The above example of the proposed embedding- 
reduction method demonstrates the advantage of ana- 
lytical solutions. Not only could we find some solutions 
of the specific type (|3"3"|) but we actually were able to ex- 
haust all solutions of this type. This also implies that 
for parametrization of flow velocity field, like e.g. in the 
blast wave model for RHIC fireball, there are only very 
limited choices for the flow profile ansatz. 
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V. SUMMARY AND DISCUSSION 

In summary, a general framework for the analytical 
treatment of RIHD equations has been developed. The 
method features a separation of longitudinal and trans- 
verse expansions, as inspired by RHIC phenomenology 
After the separation, the longitudinal and transverse ra- 
dial flows are embedded utilizing lower-dimensional solu- 
tions. The remaining equations are found to be reducible 
to a single constraint equation for transverse radial flow 
velocity field only, which can be solved completely for a 
certain ansatz for the velocity field. All solutions with 
power-law dependence on proper time and transverse ra- 
dius have been found. 

We now discuss various possible extensions of the 
present approach. 

Nontrivial longitudinal embedding: In the current work 
the longitudinal flow is embedded with the Hwa-Bjorken 
solution. It would be very interesting to try embed ding 
the newly found (1+1)D solutions in e.g.[23l [24l l25l [26f 
with more realistic longitudinal expansion for RHIC 
which would be useful for studying elliptic flow in the 
forward/backward rapidity and their correlation[33l|. 

Solutions with non-power-law transverse expansion: It 
would also be interesting to a test more nontrivial ansatz 
for the embedded transverse flow. For example we know 
from numerical calculations of radial flow [34] in central 
collisions at RHIC that the radial velocity field may be 
parameterized as v p w /(t)t"/t with f(r — ► 0) — > and 
/(t >> 1) — > 1. Such parametrization can be cast into 
the derived velocity equation (|28p to find possible solu- 
tions. 

Small deformation and elliptic flow: The analytic 
treatment of transverse elliptic flow is difficult. One ap- 
proximate method may be to introduce a parametrically 
small deformation of the matter field (with a certain ec- 
centricity parameter e 2 ) on top of an exact solution with 
transverse radial flow and using linearized hydrodynam- 
ics equations to investigate possible universal relations 
between the finall y de veloped velocity field anisotropy v 2 
and the initial e 2 |35j . 

Transverse elliptic flow embedding: Another possibil- 
ity to seek exact solutions with transverse elliptic flow 
is to use instead of (p, <fi) certain hyperbolic coordinates 
which by definition incorporate elliptic anisotropy, see 
Appendix of [361 ] for an example of such curved coordi- 
nates which may be used to develop a similar embcdding- 
rcduction procedure describing transverse elliptic flow. 
Another possibility will be combining certain conformal 
transformations with hydrodynamics equations to de- 
grade the elliptic geometry back to a spherical one. 

2D Hubble embedding: In all the previously discussed, 
we have chosen to embed (1+1)D Hubble flow for the 
(i, z) — ► (t, rj) part, due to an emphasis on RHIC evolu- 
tion. Theoretically, one can also embed a (1+2)D Hubble 
flow for the (t, x, y) — > (i, p, 4>) — ► (r p , rj p , (j>) part, and can 
eventually reduce the equations to a velocity equation 
with two variables (r p , z) in exactly the same manner as 



before. 

Acknowledgments 

The work is supported by the Director, Office of 
Energy Research, Office of High Energy and Nuclear 
Physics, Divisions of Nuclear Physics, of the U.S. De- 
partment of Energy under Contract No. DE-AC02- 
05CHII23I. JL is grateful to Edward Shuryak for very 
helpful discussions. 

Appendix A 

In this Appendix we give the detailed derivations lead- 
ing to the solutions in Section IV with the transverse 
velocity ansatz (|33l) . 

We first evaluate the derivatives ^ and ^ with J 7 , Q 
given in (|3"4"|) . Again we will make use of v p . p = 
v p - C/p and v p . t = v p -B/t for the velocity ansatz (|33[) . 

The result for |£ is 

dp 

v p ( _ f _ 2 . 

-d~p = ,(l-,)p^(l-v^ X V lPTV " {1 - V ^ 
+/ 2 r 2 [(1 + C)^ + (C-1)](1-^) 

+hrpv p + / 4 t 2 [(3C + l)v 2 p + (C- 1)] J 

= , ( i-,)Ai-ffl x {w +i)1 ^ 

+ l-h] Prv 3 p + [2/ 2 + / 4 (3C + 1)] r 2 v 2 p 

+ [h + h] prvp + [(/ 2 + U)(C - I)] r 2 | (Al) 

with coefficients /i,2,3,4 given by 

/1 = -[B + (1 - AB)v + 2Bv 2 ] x (2C) 

f 2 = -[C+(l-C)u] 

h = Bx {AC) x (I - 2v) 

U = C X (1 - 2v) (A2) 

The result for |£ is 

Or 

9Q Vp \ - -2n 

97 = Ki-^(i-^ x r^ ( p) 

+g 2P 2 [(I + B)v 2 + (B - I)] (1 - v 2 p ) 
+93Tpv p + g 4 p 2 [(3B + l)v 2 p + (B- I)] J 

+ [- gi ] prv 3 p + [2g 2 + g 4 (3B + I)] p 2 v 2 p 

+ [.9i + .93] prvp + [(ff 2 + 94)(B - I)] p 2 ) (A3) 
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with coefficients 31,2,3,4 given by 

gi = -[-C-v + 2Cv 2 \ x (2B) 
g 2 = -[(B-l)u] 
.93 = (4B) x C x {2v - 1) 
.94 = £ x (2i/- 1) 



(A4) 



Now combining the results into (|28|) we obtain the fol- 



lowing (with w p = At 

dF_dG_ 

dp dr 

0=X[t, P ] -- 



p c substituted in) 



f l) J 4 4 ]r 4S p 4C + 2 

/i)A 3 ]r 3S+l p 3C+l 



1/(1 — v)p 2 T 2 (l — V 2 

[-MC + 1)A 4 ]t 

+[(91 - 

+ [(2/ 2 + / 4 (3C+l))A 2 ]r 
+ [(-2g 2 -g 4 (3B 

+[(h + f3-gi-g 3 )A}T B 

+ [(/ 2 + /4)(C-l)]r 2 
+ h(.92 + .94)(S-l)]p 2 



J[r,p]=0 



23+2^20 



l))A 2 ]r 2B p 2C+2 



-i p c+i 



(A5) 



Again one can test the correctness of the above equation 
by using the Hwa-Bjorken (A — 0) and the 3D Hubble 
(A = 1, B = -1, C = 1) solutions. 

In equation (|A5[) . terms with various powers of r, p 
(and only power terms) appear in X[t,p]: to make all of 
them, either mutually cancel (among terms with exactly 
the same r, p powers) or vanish by respective coefficients, 
to eventually zero is quite nontrivial. A thorough sorting 
of the sequences of r, p powers can exhaust all possibili- 
ties to satisfy the algebraic equation (|A5|) . 

To see how this actually works, we give one concrete 
example. Let's consider the case when B > and C^O: 
this implies that for the exponents of r we have AB + 
2 > AB > 2B, AB + 2 > 3B + 1 > B + 1, AB + 2 > 
2B + 2 > 2B, 2B + 2 > 2 and 2B + 2 > B + 1. So the 
term [-/ 2 (C + 1)A 4 }t 4B+2 P 4C in Eq.JMJ can NOT be 
cancelled by any other one and has to vanish by itself: 
this leads to 



f 2 (C + l)A 4 = Q 



(A6) 



which in turn gives three possibilities f 2 — or C = — 1 
or A = 0. In this example we follow C = — 1 (other 
choices lead to other solutions). With C = — 1 we notice 
again that the term [— (g 2 + gi){B - l)]p 2 can NOT by 
cancelled by any other remaining terms and thus shall 
vanish by itself: this leads to 



( 52 +.94) (5-1) = 



(A7) 



which again has two possibilities B = 1 or g 2 + g± = 0. 
Now we choose to follow B = 1: with this choice the 
remaining terms are significantly simplified and finally 
lead to two equations about the coefficients: 

2g 2 A 4 + ( 5l - /i)A 3 + 2(/ 2 - f 4 )A 2 = 
-2(02 + 2 9i )A 2 + {h +f 3 - 9l - g 3 )A - 2(/ 2 + / 4 ) = 

(A8) 



It can then be verified that the only solution is A = 1 for 
arbitrary v. 

Of course there are many but finite number of combi- 
nations that one can follow to check one by one. Note 
not all possibilities appearing initially can finally lead to 
a solution: there are only four variables A,B,C,v and in 
most cases it turns out contradiction occurs at the end 
which means no solution. After a tedious examination 
we have found all possible solutions as listed in Section 
IV, and there is no more solution of the power law ansatz 
type as in (|33j). 



Appendix B 

In this Appendix we use (1+1)D ideal relativistic hy- 
drodynamics to demonstrate the reduction method in a 
more explicit manner. The hydrodynamics equations are 
(in (t,z) coordinates) 



[d t + vd. 



= -dzv-^dt+vd^-) (Bl) 



d z p vd t p 



(B2) 



e + p e + p 

In the above v is the spatial velocity dz/dt and j v = 
l/Vl — v 2 . The energy density e and pressure p shall be 
related by the E.o.S, which we use in a slightly different 
way. We introduce the enthalpy density w = e + p and 
the speed of sound c s = ^Jdp/de (which can be deduced 
from E.o.S), and use the following relations 



1 c z 

de = -dw , dp — ^-rrdw 

1 + c 2 1 + c 2 

to re-write the hydrodynamics equations into: 



(B3) 



<9t[ln(u;)] + 1; 9 2 [ln(w)] 
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[d x v + vd t v] (B4) 



vd t [\n(w)]+d z [\n(w)] = -^-[vd x v + d t v] (B5) 

with £ = c 2 /(l + c 2 ). From these two equations we can 
obtain <9 t [ln(u>)] and i9 2 [ln(u;)]: 

d t [ln(w)} = X[v, d t v, 3 z v] = X[t, z] 



-it 



- (1 - £)v 2 ]d z v + [2£-i\vd t v 



(B6) 



d z [ln(w)]=y[v ) d t v,d z v]=y[t > z] 



[(l-0-£v 2 ]d t v+ [l-2£]vd a 



(B7) 



The necessary and sufficient conditions for the above set 
of equations to be soluble is the following: 



d z x[t,z]-d t y[t,z] = 



(B8) 
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Thus we have reduced the original hydrodynamics equa- 
tions into a single but second order differential equation 
for the velocity field v(t, z) only. With the above satis- 
fied, the matter field is given by 

/, x f* dt' X[t' ,z]+ P dz'y[t,z'] (d m 

w{t, z) = w ■ e Jt o z o (B9) 

with wo its value at arbitrary reference point (to, Zq). 

In the case of a linear E.o.S as the one in ([3]), the speed 
of sound c s and thus £ are constants independent of e or 
p, and we can further simplify the reduced equation (|B8|) 
for velocity field into the following: 

[(l-Z)-Zv 2 ](l-v 2 )(d t d t v) 
+ [(2 - 30 - Zv 2 }(2v)(dtv) 2 

~m-l) + ^-l)v 2 ](2v)(d z v) 2 

+2(1- 2^(1 ~v 2 )v(d t d z v) 

+2(1 - 2^(1 + 3v 2 ){d t v)(d z v) = (BIO) 

A similar scheme can be carried out for curved coor- 
dinates like (t, 77) in a straightforward way. We no- 
tice a similar imp lementation using light-cone variables 
z± = t ± z in [24j . 

We emphasize that while the above procedure seems 
somewhat trivial in (1+1)D, its realization is much more 
nontrivial and involved in (1+3)D. We also point out that 
the reduced equation (|B8p for the velocity field (or the 
simplified one in case of linear E.o.S) has to be satisfied 
by all solutions to the (1+1)D hydrodynamics equations. 
In Appendix C we give a nontrivial and involved exam- 
ple from the recently found Nagy-Csorgo-Csanad (NCC) 
solutions [231] (which also include (1+1)D Hwa-Bjroken as 
a special case) to show the correctness and usefulness of 
the derived velocity equation. 



Appendix C 

In this Appendix we show that the NCC family of ana- 
lytic solutions in [23[ for 1-D ideal hydrodynamics equa- 
tions with a linear E.o.S can also be deduced by subject- 
ing their velocity field ansatz to the reduced equations 
(|B10[) we derived. With the resulting velocity field we 
also show the matter field of NCC solutions is indeed 
given by (|B9j) . 

The velocity field ansatz of NCC solutions is the fol- 
lowing (for inside-light-cone region, i.e. \z\ < \t\): 

(t + z) x - {t-z) x 



tanh[Ary] = 



(t + z) x + (t- z y 



v : 



2 l t-z l 



(CI) 

with A some constant. With the above, we obtain the 
following relations for the derivatives: 

d t v = A(l - v 2 )d t r, , d z v = A(l - v 2 )d zV 

d t d t v = A(i - v 2 )[{d t d tr] ) - 2\v{d t i 1 ) 2 ] 

d z d z v = A(l - v 2 )[(d z d zV ) - 2Xv(d z r 1 ) 2 } 

d t 8 z v = A(l - v 2 )[(8 t d zV ) - 2Xv(d t r 1 )(d zV )] (C2) 



Substituting the above derivatives into our velocity equa- 
tion ([B1Q]) . we obtain the following: 

= A(1-tj 2 ) 2 x 

[(1 - - tAidtdtr!) + [-i + (1 - 0A{d z d z r,) 
+2(1 - 2ti)v(d t d zV ) + 2(1 - 2£)A(1 + v 2 )(d t ri)(d z ri) 
+2(l-20Xv[(d t r 1 ) 2 + (d z r 1 ) 2 }) (C3) 



After evaluating the derivatives of 77 in the above, we 
obtain: 



2A(1- v 2 



(l-20-(l-A)-[(teK -it +z 2 )v + {tz)] = 



(C4) 



{t 2 - z 2 ) 2 
We find three classes of solutions 

• £ = 1/2 with arbitrary A; 



A = 1 with arbitrary £ (which is nothing but the 
Hwa-Bjorken solution); 



• \{tz)v 2 - (t 2 + z 2 )v + (tz)] = which yields only 
one meaningful solution with v = z/t, but this is 
just the A = 1 solution. 

These cover all (1+1)D NCC solutions found in 
Note their parameter k from E.o.S e = np is related to 
our E.o.S parameter £ = c 2 /(l + c 2 ) by k = (1 — 0/0 

Next we examine the matter field corresponding to the 
velocity field solutions: 

• for £ = 1/2 case: we have X = (-\)[2t/(t 2 - z 2 )] 
and y = (— A)[— 2z/(t 2 — z 2 )], which via our 

equation (|B9|) gives w = w [ ^~^ ] A ; 



• for A = 1 case: we have v — z/t and thus 
X = ^{2t/ (t 2 - z 2 )] and y = ^[-2z/{t 2 - 
z 2 )], which via our equation (|B9p gives w = 

-tl-z 2 a i 1/(2-2?) 



The two cases can be combined into a single form: 



t 2 



which is the same as obtained in 12311. 



(C5) 
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